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ABSTRACT 

Pre-transitional disks are protoplanetary disks with a gapped disk structure, potentially indicating 
the presence of young planets in these systems. In order to explore the structure of these objects 
and their gap-opening mechanism, we observed the pre-transitional disk V1247 Orionis using the Very 
Large Telescope Interferometer, the Keck Interferometer, Keck-II, Gemini South, and IRTF. This 
allows us spatially resolve the AU-scale disk structure from near- to mid-infrared wavelengths (1.5 to 
13 /im), tracing material at different temperatures and over a wide range of stellocentric radii. Our 
observations reveal a narrow, optically-thick inner-disk component (located at 0.18 AU from the star) 
that is separated from the optically thick outer disk (radii > 46 AU), providing unambiguous evidence 
for the existence of a gap in this pre-transitional disk. Surprisingly, we find that the gap region is 
filled with significant amounts of optically thin material with a carbon-dominated dust mineralogy. 
The presence of this optically thin gap material cannot be deduced solely from the spectral energy 
distribution, yet it is the dominant contributor at mid-infrared wavelengths. Furthermore, using 
Keck/NIRC2 aperture masking observations in the H, K', and V band, we detect asymmetries 
in the brightness distribution on scales of ~ 15-40 AU, i.e. within the gap region. The detected 
asymmetries are highly significant, yet their amplitude and direction changes with wavelength, which 
is not consistent with a companion interpretation but indicates an inhomogeneous distribution of the 
gap material. We interpret this as strong evidence for the presence of complex density structures, 
possibly reflecting the dynamical interaction of the disk material with sub-stellar mass bodies that 
are responsible for the gap clearing. 

Subject headings: stars: pre-main-sequence — planetary systems: protoplanetary disks — accretion, 
accretion disks - techniques: interferometric 
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1. INTRODUCTION 

Planets are believed to form in the circumstellar 
disks around you ng stars, either th rough a process 
of core accretion ([Pollack et al.l 11996ft or gravita tional 
insta bilities in the more extended disk regions (|Bossl 
2000). A particularly interesting phase in this pro- 
cess starts when the newly-formed planetary bodies 
have gained sufficient mass to interact with the am- 
bient disk mat erial and affect the disk struct ure sig- 
nificantly (e.g. iPaardekooper fc Mellemal I2004P) . Po- 
tential candidates for disks that might have been dy- 
namically aff ected by planeta ry bodies are the transi- 
tional disks (S trom et al.lll989T ) and pre-transitional disks 
(Esp aillat et al.l I2007T ). These objects exhibit a strong 

1 Based on observations made with the Keck observatory 
(NASA program ID N121IV and NOAO program ID N121N2), 
Gemini South (NOAO program ID GS-2011B-Q-19), and with 
ESO telescopes at the Paranal Observatory (ESO program IDs 
088.C-0868(A) and 088.C-0763(A+B)). 



far-infrared (> 20 fim) excess, but have a significantly 
reduced near-infrared (NIR) to mid-infrared (MIR) ex- 
cess compared to classical T Tauri disks. This reduced 
excess emission indicates that the inner-most disk re- 
gions contain only optically thin gas and dust (tran- 
sitional disks) or exhibit an extended gap, which sep- 
arates the optically thick inner disk from the outer 
disk (pre-transitional disks). The inner "holes" and 
gaps observed in transitional/pre-transitio nal disks could 
be caused by d isk-planet interaction (e.g. lOuillen et al.l 
l2004tlK"raus fc IreTa nd 2012), but several alternative disk 
clearing sce narios have been proposed, including grain 
growth (e.g.lDullemond & Dominik 2005 ), magne torota- 
tional instabili ties dChiang fc Murrav-Clavl[20 07\ pho- 
toevaporation ((Alexander fc Armitagell2007D. and trun- 
catio n by close-in stellar companions ( Ireland fc Krausl 
(2008ft . 

Most of the aforementioned processes take place in the 
inner few astronomical units (AU) around the central 
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star, corresponding to angular scales p < 0.01 ", even 
for the nearest young stars. Given that these scales are 
not accessible with conventional imaging techniques, ear- 
lier studies relied mostly on the modelling of spatially 
unresolved constraints, in particular the spectral energy 
distribution (SED). However, these modelling techniques 
suffer from well-known degeneracies - grain temperature, 
for instance, is sensitive to both its distance from the illu- 
minating star, the pa r ticle size, and the dus t composition 
(jThamm et all 119941 : iVinkovic et"aH I2003Q . Therefore, 
we initiated an observational campaign to study these ob- 
jects using a multi-wavelength interferometry approach, 
where we combine near- and mid-infrared interferometry 
constraints. The data set includes long-baseline interfer- 
ometric data from the Keck Interferometer (KI) and Very 
Large Telescope Interferometer (VLTI), single-dish inter- 
ferometric data from Keck & Gemini South, and spec- 
troscopic data from IRTF and the ESO 3.6 m telescope. 
Given the wide wavelength coverage (1.5 — 13/im), our 
interferometric observations probe a wide range of dust 
temperatures and material located over a wide range of 
stellocentric radii, from sub-AU to tens of AU. 

Here, we report first results from our campaign, fo- 
cusing on the pre-transitional disk object V1247 Orionis 
(HD290764, PDS192). V1247 0ri is likely a member 
of the Orion OB lb association (|Schild fe Cowlevl 119711: 
iGuetterl I1981D and the Alnilam cluster (also known as 
"e Orionis cluster" and "Col linder 70"), whose age ha s 
been estimated to 5-10Myr (jCaballero fe SolanoH 2008). 
There is also considerable debate in the literature con- 
cerning the precise spectral classific ation of V12470ri , 
ranging fr om spectral types of F OV ()Vieira et al.ll2003T) 
to A5III (|Schild fe Cowlevl fl97ll) . which has caused us 
to re-eval u ate th e spectral classification in our study. 
Caballcro (2010) reported the observation of two iso- 
lated, but deep UXOri-like occultation events, which 
were attributed to occultation from disk material in the 
inner disk regions. The SED of V12470ri shows a signif- 
icantly red uced exce ss emission in the wavelength range 
3 — 15 (im (Caballcro 2010), revealing similarities with 
other pre-transitional disk objects. However, V1247 Ori 
is of considerably earlier spectral type than most other 
pre-transitional disks, which makes this system an in- 
teresting laborary to study disk evolution and planet- 
formation in the intermediate-mass regime (~ 2Mq). 

In the following, we present our extensive observa- 
tional data set (Sect. [2]) and our results, including a 
re-evaluation of the spectral classification of the cen- 
tral star (Sect. I3.1|) . the SED and line diagnostics 
(Sect. 13. 2[) . our interferometric constraints on the disk 
structure (Sect. I3.3[) . and a global model based on our 
multi-wavelength interferometry and SED constraints 
(Sect. I3.4J) . Finally, we discuss the implications of our 
work on the understanding of the pre-transitional disk 
phase (Sect. ID) and summarize our conclusions (Sect. [5]). 

2. OBSERVATIONS 

2.1. IRTF/SpeX+BASS near- and mid-infrared 
spectroscopy and archival spectro-photometry 

An important first step is to build the SED of 
V1247 0ri in the relevant near- to mid-infrared regime. 
For this purpose, we emplo yed the SpeX spectrograph 
(l-5^m: lRavner et al.ll2003h and The Aerospace Corpo- 



ration's Broad-band Array Spectrograph System (BASS; 
3-14 /mi), which are mounted at NASA's Infrared Tele- 
scope Facility (IRTF). The observing dates and observ- 
ing modes are summarized in Table [TJ together with our 
other observations. 

The SpeX spectra were recorded using the echelle grat- 
ing in both short-wavelength mode (SXD, 0.8-2.4 /im) 
and long-wavelength mode (LXD, 2.3-5.4 /urn) using a 
0.8" slit. The spectra were corrected for telluric ex- 
tinction and flux calibrated against the calibrator star 
HP 34203 using the Spextool data re duction package 
(jVacca et all 120031: iCushing et all 120041 ). In addition to 
the 0.8"-slit spectra, we also recorded data with a wide 
3" slit, which allows us to retrieve the absolute flux lev- 
els. VI 247 Ori was observed six times in total with SpeX 
(see Table [1]), allowing us to detect possible mid- infrared 
variability. 

Our BASS short-wavelength (2.9-6 /im) and long- 
wavelength (6-13.5 /im) observations cover the full L, M, 
and N bands with a spectral resolution ranging from 
about 30 to 125. For extracting the BASS data, we used 
the Spextool data reduction package. The calibrator star 
was HD 29139. 

The UBVRI photome try was taken from the compila- 
tion by ICaballerol (l2010f). which co n tains measurements 
from iSharplessI (Il962h, iFuiii et al.1 (12002D. the I nfrared 
Astronomical Satellite (IRAS. iHelou fe Walker! r 1 988). 
the Midcourse Space Exverimen t (MSX. | Egan et al] 
|2003[ L and the Akari observatory (jlshihara et al.ll2010f ). 
In addit ion, we added JHK band photometry from 
2MASS (ISkrutskie et all l200l . Besides the ASAS V- 
band data published in ICaballerol (|2010l) . no new ASAS 
observations are available. The Spitzer IRS data were 
obtained from the Spitzer heritage Archive Post-Basic 
Calibrated Data (PBCD). All spectra were low resolu- 
tion spectra, with the deviant end points of each spec- 
tral segment trimmed off. The resulting SED is shown 
in Fig. n 

2.2. KI near-infrared interferometry 

The Keck Interferometer allowed us to resolve the dis- 
tribution of the hot (T > 1000 K) circumstellar material 
around V12470ri, forming an effective projected base- 
line length B p of 70.1 m towards position angle (PA, mea- 
sured East of North) 27°. The observation was conducted 
on 2011 November 5 and used the A STRA beam com- 
bination instrument (jPott et al.l 2009) in its sclf-phase- 
referencing mode (V2-SPR), covering the full K band 
with spectral resolution R = A/AA = 2000. The science 
target observations were accompanied by observations of 
the calibrator stars HD 41794 and HD 37331, which al- 
lowed us to correct for atmospheric and instrumental 
effects. For the calibrator sizes we adopt uniform disk 
(UD) diameters of 0.213 ±0.015 milliarcsecond (mas) for 
HD 41794 (computed with ASPRO) and 0.144 ± 0.1 mas 
for HD 37331 (computed with ASPRO). 

From these observations, we derived a spectrum as 
well as wavelength-differential visibilities and phases. 
However, neither the spectrum nor the interferomet- 
ric observables show any significant spectral features or 
wavelength-differential features. Therefore, we use only 
the recorded broadband visibility information for the 
modeling described in the following sections. 



Resolving the gap and AU-scale asymmetries in the pre-transitional disk of VI 247 Orionis 



TABLE 1 
Observation log 



Instrument 


UT Date 
[YYYY/MM/DD] 


Mode 


Seeing 
[arcsec] 


Calibrator 


IRTF/SpeX 


2010/02/28 


SXD 
LXD 
Prism 


0.6 


HD 34203 
HD 34203 
HD 34203 


IRTF/SpeX a 


2011/03/04 


SXD 


0.8 


HD 34203 


IRTF/BASS b 


2011/10/17 




0.7 


HD 29139 


IRTF/SpeX c 


2011/10/18 


SXD 
LXD 


0.6 


HD 34203 
HD 34203 


KI/ ASTRA 


2011/11/05 


V2-SPR 


0.6 


HD 41794, HD 37331 


Gemini S/T-ReCS 


2011/11/23 


Si-2 
Si-3 
Si-5 


0.5 


HD 36167 
HD 36167 
HD 36167 


VLTI/AMBER 


2011/12/09 


UR-HK 


1.1 


HD 37331, HD 36059 


VLTI/MIDI 


2011/12/11 


Prism 


0.9 


HD 39400 


VLTI/MIDI 


2011/12/12 


Prism 


0.7 


HD 39400 


Keck-II/NIRC2 


2012/01/08 


K' 


0.1 


HD 37634, HD 38406 


Keck-II/NIRC2 


2012/01/10 


H 


1.1 


HD 37331 






V 


1.3 


HD 37634, HD 38406 


IRTF/SpeX 


2012/02/26 


SXD 
LXD 
Prism 


0.9 


HD 34203 
HD 34203 
HD 34203 


IRTF/SpeX 


2012/09/12 


SXD 
LXD 
Prism 


0.8 


HD 34203 
HD 34203 
HD 34203 


IRTF/SpeX 


2012/11/04 


SXD 
LXD 
Prism 


0.4 


HD 34203 
HD 34203 
HD 34203 



a SXD only due to clouds (1-mag variations in sky transparency). 

b Exceptionally low precipitable water vapor allowed observations through the 6 (im 

water vapor band. 

c No Prism observations due to thin cirrus. 
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Fig. 1.— SED of V1247 0ri, compiled from our IRTF/SpeX (epoch 2010 February 28) and IRTF/BASS (epoch 2011 January 17) 
observations and data from the literature. The horizontal error bars indicate the ban dpass of the e mployed photometric filters. The 
photospheric emission (blue and magenta curve) is represented with the SYNTHE model (Kurucz 1993) for a F0V star with T e g = 7250 K 
and log g = 4.5. 



2.3. VLTI/AMBER near-infrared interferometry 

Additional NIR long-baseline interferometric observa- 
tions were obtain ed on 2011 Decem ber 9 with the AM- 
BER instrument ()Petrov et al. 2007), improving the uv- 
plane coverage of the KI interferometry. Using the VLTI 
8.2 m unit telescope (UT) triplet UT2-UT3-UT4, the ob- 
servations sample projected baseline lengths (position 
angles) of 43.4m (221°), 62.4m (289°), and 88.3m (82°). 
Using optical fibers, AMBER provides spatial filtering, 
resulting in an effective field-of-view of 60 mas. For the 



beam combination we used AMBER's LR-i/i^-modc, 
which provides a spectral dispersion of R = 35 and 
covers wavelengths between 1.5 and 2.5 fim {H- and K 
band). We recorded 5000 interferograms with an inte- 
gration time of 26 ms for V12470ri and the calibrators 
HD 37331 and HD 36059 (UD 0.549 ± 0.038 mas, AS- 
PRO). Data from another calibrator (HD 36811) were 
rejected, since this object was found to be a close binary 
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Fig. 2. — Squared visibility amplitudes (top) and closure phases 
(bottom) in the H- and K band measured with VLTI/ AMBER. 
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Fig. 3. — Comparison of our VLTI/MIDI spectrum with the 
mid-infrared spectrum recorded with the IRTF/BASS. The shaded 
area marks the location of the atmospheric ozone absorption bands, 
which might introduce additional scatter in the spectrum. 

systeixQ. 

From the data we extract wavelength-dependent visi- 
bilities and closure phases using the a mdlib software (Re- 
lease 3; EtuIireEaD[2002[ChSreEaD[lQp). In order to 
minimize the effect of residual telescope jitter, we follow 
the standard AMBER data reduction procedure and se- 
lect the interferograms with the 10% best signal-to-noise 
ratio. In addition, we reject scans where the optical path 
delay exceeds 4 /im. The wavelength calibration was done 
using the telluric absorption ba nds between the J/ H and 
H/K band (see Appendix A in lKraus et"ai1 l20"09b). The 
final calibrated observables are shown in Fig. [3J 

2 Given that our companion detection for HD 36811 might be 
of interest for science programs on Am-type stars, we determined 
the binary parameters and obtained a separation p = 2.0 ±0.2 mas 
with position angle PA=80 ± 2° (epoch 2011-12-09) and flux ratio 
F B /F A = 0.44 ± 0.01 in the H/K band. 
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Fig. 4. — Visibilities derived from our VLTI/MIDI observations. 
The shaded area marks the location of the atmospheric ozone ab- 
sorption bands, which might introduce additional scatter in the 
spectrum. 



2.4. VLTI/MIDI mid-infrared interferometry 

Our mid-infrared long-baseline interferometric obser- 
vations were obtained w ith the VLTI/MIDI instrument 
(jPrzvgodda et al1 l2003). The employed telescope base- 
lines were UT2-UT3 (2011 December 11, B p = 46.6 m, 
PA=46°) and UT1-UT4 (2011 December 12, B p = 
125.0m, PA=63°), i.e. probing similar position angles, 
but significantly different spatial frequencies. The obser- 
vations on V1247 0ri were bracketed with observations 
of the calibrator star HD 39400 (UD 2.387 ±0.138 mas, 
computed using CalVin). 

In order to extract the mid-infrared correlated 
fluxes, we employ the MIA+EWS-assisted (| Jaffei 120041: 
iLeinert et al. 2004) fai nt-source red uction procedures de- 
scribed in Kishim oto et al.1 (J2011I ). Together with the 
MIDI interferograms, we also recorded photometry files 
on the science star and calibrator, which allowed us 
to retrieve the absolute flux photometry. Comparing 
the MIDI total flux spectrum with a BASS spectrum 
recorded two months earlier (Fig. [3]) reveals that the 
spectra are fully consistent within the measurement un- 
certainties and we find no indications for mid-infrared 
flux variability. Therefore, we decided to use the BASS 
spectrum in order to convert the measured correlated 
fluxes to visibility amplitudes (Fig.U]). 

The wavelength range covered by MIDI contains also 
two hydrocarbon features (Sect. I3.2J) . However, we de- 
tect no significant change in the visibility at the wave- 
length of these organic bands (Fig. 2]), suggesting that the 
hydrocarbon grains are located on similar spatial scales 
as the optically thin dust grains. However, observations 
with shorter baselines and with higher signal-to-noise ra- 
tio will be required for a detailed quantiative analysis. 

2.5. Gemini/T-ReCS mid-infrared speckle 
interferometry 

Our MIDI observations show that V12470ri is al- 
ready strongly resolved on the shortest VLTI UT baseline 
(B p = 46.6 m), indicating the presence of an extended 
emission component that is not properly constrained by 
our long-baseline interferometric observations. There- 
fore, we conducted complementary observations using 
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Fig. 5. — 2-D power spectra, derived from our Gemini/T-ReCS 
speckle observations using wavelength bins around 8.74 /im (Si-2, 
top), 9.69 /im (Si-3, middle), and 12.66 (im (Si-5, bottom). For each 
filter we recorded two calibrated data sets on V1247 Ori (left and 
right), which confirm the detected object elongation independently. 
The contours mark visibility levels of 0.7, 0.8, and 0.9 in the derived 
power spectra. North is up and East is left. 



the T-ReCS mid-infrare d imager (jTelesco et all 119981 : 
IDe Buizer fc Fisher 20Q5|) that is mounted on the Gem- 
ini/South 8.2 m telescope. Employing short integration 
times of 0.2 s and a speckle interferometry analysis ap- 
proach allow us to effectively freeze the atmospheric per- 
turbations and to extract interferometric visibilities and 
phases for effective baseline lengths B p < 5.5 m that 
can be directly combined with our long-baseline inter- 
ferometric data. The data were recorded on UT 2012 
November 23 under exceptional atmospheric conditions 
using the T-ReCS narrowband filters Si-2 (A c = 8.74 /im, 
AA = 0.39 /im, where A c denotes the central wave- 
length of the filter and A A the spectral bandwidth), Si-3 
(A c = 9.69 /im, AA = 0.46 /im), and Si-5 (A c = 11.66 /im, 
AA = 0.56 /im). Compared to the other two channels, 
the Si-3 data exhibit a reduced signal-to-noise ratio, re- 
flecting the low flux of the object at this wavelength and 
the presence of the atmospheric ozone absorption bands. 
As a calibrator star, we employed HD 36167. 

The power spectrum a nalysis method is described in 
more detail in Eisner et al. (2009). Speckle observations 
in the mid-infrared are rather sensitive to variations in 
the thermal background, which can potentially induce a 
bias in the absolute visibility level. In order to correct 
this problem, we renormalize the measured visibility pro- 
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Fig. 6. — Squared visibility amplitudes derived from our 
Keck/NIRC2 aperture masking observations. 

files on the shortest spatial frequencies. 

The derived 2-D visibilities show consistently an ob- 
ject elongation (Fig. [5j that we will further quantify in 
Sect. I3.3.1I using model-fitting techniques. The validity 
of this elongation is also confirmed by the fact that the 
T-ReCS data were recorded with deactivated field ro- 
tator, which would introduce a significant PA rotation 
in the derived 2-D visibilities for any telescope-induced 
artefacts. 

2.6. Keck-II/NIRC2 near-infrared aperture masking 

Using spatial filtering by single-mode fibers, our AM- 
BER observations are insensitive to structures located 
at separations p larger than ~ 30 mas from the central 
star. In order to probe larger spatial scales, we em- 
ployed aperture masking observations with the Keck-II 
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Fig. 7. — Top: Our Keck/NIRC2 .ff-band (left column), /{"'-band (middle column), and L'-band (right column) aperture masking 
observations reveal strong non-zero phases that we fit with a companion model (solid line). The 2-D Fourier phases have been projected 
perpendicular to the position angle of the best-fit companion model. Bottom: x 2 -maps derived by fitting a companion model. The best-fit 
companion positions do not agree for the different wavele ngth bands, indicating that the measured asymmetries are not caused by a simple 
point source, but complex disk asymmetries (Sect. [3T373J) . 

10 m telescope and the NIRC2 instrument. Using a non- 
redundant 9-hole pupil mask allows us to derive high- 
precision interferometric phases that are well-suited to 
search for faint companions. The observations were con- 
ducted on 2012 January 8 and 10 with an integration 
time of 1 s and with the adaptive optics loop closed. The 
observations on 2012 January 8 were conducted under 
excellent atmospheric conditions using a X'-band (1.96- 
2.29 //m) filter. We executed an observing sequence of 
four target-calibrator pairs with a total on-source in- 
tegration time of 960 s. On January 10, four addi- 
tional sequences were recorded using an H-band (1.5- 
1.8/im), L'-band (3.43-4.13^m), and again a -ftT'-band 
filter. The observed calibrator stars were HD 37331 (H 
band), HD 37634 (K'/L 1 band), and HD 38406 (K'/ L'- 
band) . 

The NIRC2 data were reduced using a data reduction 
pipe line that was used already in various earlier s tudies 
(e.g. Ilreland fc Krausll200cl iKraus fc Irelandll2"o"l2T ) . pro- 
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viding absolute calibrated visibilities (Fig. [5]) and closure 
phases (Fig. [3 top). 

2.7. Archival HARPS visual spectroscopy 

Given the considerable uncertainties on the spectral 
type of V1247 Ori, we obtained HARPS high-resolution 
spectra from the ESO archive. The 37 exposures were 
reduced with the HARPS DRS pipeline v3.0 and pro- 
vide a resolution R = 120, 000. The data were recorded 
between 2008-11-11 and 2008-11-13 in the course of a 
radial velocity survey (ESO program ID 082.D-0833(A), 
PI: Boehm). Based on this original science objective, the 
data were recorded without associated telluric calibra- 
tors. Therefore, we were not able to remove atmospheric 
features and thus focused our analysis on the spectral re- 
gions around 380-460 nm and around the Ha line, which 



Wavelength [nm] 

Fig. 8.— Comparison of the HARPS spectrum of V1247 Ori (3rd 
from top) for the spectral window between 380 and 460 nm with 
synthetic spectra (Munari ct al. 2005) for effective temperatures 
of 10,000 K (top), 8500 K (2nd from top), 7250 K (4th from top), 
and 6750 K (5th from top). In order to enable an easy comparison, 
we overplot the synthetic spectra (red curves) with the measured 
HARPS spectrum (blue curves). 

suffers only from very weak telluric contamination. The 
individual spectra were averaged and the continuum-level 
was then normalized by fitting a higher-order polyno- 
mial. 

3. RESULTS 

3.1. Spectral classification 

In the literature, an unusually wide range of spectral 
types has been assi gned to VI 2 47 Ori, ranging from FOV 
(|Vieira et al.1 12003) to A5III (ISchild fc Cowlevl fl97l . 
This has motivated us to re-evaluate the spectral classifi- 
cation by fitting the HARPS spectrum with the synthetic 
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Fig. 9. — HR diagram, showing the location of V124 7Qri over- 
plotted on evolutionary tracks from Brcssan ct al. (2012). We show 
a series of isochrones from 1 Myr to lOMyr in steps of 0.2 dex 
(solid lines, from top to bottom) and evolutionary tracks from 1.5 
to 3.0 Mq in steps of 0.5 Mq (dashed lines). 



TABLE 2 
Stellar parameters of V1247 0ri 



Parameter 




Unit 


Value 


Effective temperature 


T e ff 


K 


7250 


Gravity 


log 9 


cms -2 


4.5 


Metallicity 


[M/H] 




-0.5 


Microturbulent velocity 


€ 


kms -1 


2 


Rotation velocity 


"rot 


kms -1 


60 


Enhancement parameter 


[o/Fe] 







Distance 


d 


pc 


385 dz 15 


Age 




Myr 


7.4 ±0.4 


Mass 


Ah 


Mq 


1.86 ±0.02 


Stellar radius 


R* 


Rq 


2.3 ±0.3 


Radial velocity 


Vrad 


kms 1 


±20.7 ±0.5 



Note. — For details about the individual parameters and 
how they have been derived, we refer the reader to Sect. 13.11 
The radial velocity has been converted to the heliocentric sys- 
tem. 



spectra in the library of iMunari et al.l (|2005t ). which cov- 
ers photospheres over a wide range of values for effective 
temperature T c g, gravity logg, metallicity [M/H], micro- 
turbulent velocity £, rotation velocity u rot , and [a/Fe] 
enhancement. For a detailed description of the param- 
eters and the grid properties, we refer the reader to 
IMunari et al.l (|2005t ). The range of possible rotation ve- 
locities is rather coarsely sampled in the original grid. 
Therefore, we employed the unbroadened spectra and 
convolved them with the line spread function for a limb- 
darkened rotating photosphere. Our new grid samples 
the rotation velocity with a resolution of 10 kms -1 . We 
selected the best-fit model based on the % 2 -value and 
list the best-fit parameters in Tab. [U These parameters 
correspond to a spectral type F0V (based on the table 
provide d bvlKu rucz 1993;), which confirms the classifica- 
tion bv iVieira et"aLl (|2003t) . Fig. |5] shows the V12470ri 
HARPS spectrum (blue curve) overplotted with the best- 
fit template spectrum (labeled F0V) and template spec- 
tra for stars of spectral type A0V (T c g = 9520 K, log g = 
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Fig. 10.— Comparison of the SED of V1247 Ori with the SED of 
another star with a pre-transitional disk (HD 135344B; data from 
2009 May 20). 



±4.0cms" 2 ), A5V (T eB = 8200K, logs = ±4.0cms- 2 ), 
and F2V (T cff = 6890 K, \ogg = ±4.0cms" 2 ). 

As distance for V1247 0ri we adopt 3 8 5 ± 1 5 pc, 
which is the value determined by iCaballerol (|2008l ) and 
iTerrell et al.l (|2007j ) for other members of the Alnilam 
cluster. From the distance and the apparent V-band 
magnitude (mv = 9.861 ± 0.033 mag Caballcro 2010), 
we derive the absolute magntiude (My = 1.934 ± 
0.118 mag). Usi ng the colour excess (E(B-V)=0.02, 
IVieira et al.ll2003f ). and a bolometric correction for a FOV 
star (BCy = —0. 008; assuming an ag e of 6.3 Myr and a 
mass of 1.9 Mq; iBressan et al.|[2012D . we estimate the 
bolometric luminosity to L D oi = 13.43 ± 1.5 Lq. The ra- 
dius has then been derived to R+ — 2.3 ± 0.3 Rq using 
the relation L DO i = 47ri? 2 <7T c fr, where a is the Stefan- 
Boltzmann constant. 

Using the p r e-mai n-sequence evolutionary tracks from 
IBressan etall (J2012L PARSEC release VI. 0), we estimate 
the mass of V12470ri to 1.86 ± 0.02 M©, with an age 
of 7.4 ± 0.4 Myr (Fig. EJ). The quoted statistical error 
bars have been derived assuming the aforementioned lu- 
minosity uncertainties and an uncertainty of 100 K on the 
effective temperature, but do not include systematic un- 
certainties that might be associated with the calibration 
of the stellar evolutionary models. The derived age is 
consistent with the age of the Alnilam cluster (5-10 Myr, 
ICaballero fc Solanoll2008h . 

3.2. Spectral energy distribution, line emission 
identification, and variability 

The SED of VI 247 Ori shows remarkable similari- 
ties with the SED of other stars with pr e-transitional 
disks , such as H D 135344B (=SAO 2064 62. iGradv et all 
*2009tlSitko et al.pOll or HD 169142 (jGradv et aJJl2007t 
Panic et aLll2008l : iHonda et al.ll2012f ). including a strong 
infrared excess emission with a significant flux depression 
in the ~3-15/xm mid- infrared wavelength regime com- 
pared to classical T Tauri disks (Fig. [TU]). These spec- 
tral characteristics are often interpreted as evidence for a 
gapped disk structure, where the near-infrared (< 3/im) 
emission is attributed to coming from a hot inner disk 
and the mid-infrared (> 10/im) emission to colder ma- 
terial in the outer disk. 

In the K band, no lines corresponding to the CO 
overtone (Av = 2, around 2.3 /xm) or He I (2.058 /xm) 
have been detected. However, our spectra show some 
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Fig. 11. — V1247 0ri spectrum between 2 and 18 /im, showing 
strong line emission, where the strongest lines appear to be asso- 
ciated with hydrocarbon features (grey dotted lines). 
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Fig. 12.— HARPS spectrum of V1247 0ri around the Ha-line 
(black line, original resolution), overplottcd with a template spec- 
trum for a FOV photosphere (blue line). The spectrum has not 
been corrected for telluric features. 



extremely strong emission line components around 3.3, 
6.2, 7.7, 8.6, and 11.3^m (Fig. HI]), which we identify 
as hydrocarbon-related features (sometimes incorrectly 
referred to in the literature as polycyclic aromatic hy- 
droca rbon (PAH) features, e.g. Ivan Diedenhoven et al.l 
2004). We also detect accretion-tracing spectral lines, 
including Ha (Fig. [TJ]), He I, Pa7, and Br 7 (Fig. [TBI) , 
which are superposed on the underlying photospheric 
component. The equivalent width of Ha is -0.295 nm 
and -1.006 nm after we correct for the strong underlying 
photospheric absorption component (Fig. [T2l . The Ha- 
line profile is slightly asymmetric with a weak redshifted 
absorption component, showing a similar morphology as 
predicted by magnetospheric accretion models for disks 
seens closer to edge-on than face-on ([Muzerolle et al.l 
1998). 

The Br7 line is rather weak, but shows significant vari- 
ability, which is most evident between our epoch in 2010 
and the three epochs in 2012 (Fig. IT3"1) . This variability 
in Br7 seems to be correlated with the Pa 7 and He I-line 
luminosity (Fig. [IB] top), where He I even switches from 
absorption to emission between 2010 February 28 and 
2012 February 26. Furthermore, there appears to be a 
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Fig. 13. — SpeX spectra covering the He I and Pa7 line (top 
panel) as well as the Br7 line (bottom panel) at four epochs. 




10" 



E 
5 



10" 



<t, ■ ■ ' 1 




^H^ 


a *' 55 * iL 


; o 960306 


Fuji! 2002 ^"^ 


A 981030 


2 MASS ^f I^l, 


— 080326 


spltzer H I 


— 100228 


SpeX 1ILJ a' 


• 111017 


BASS l] ThJ\ 


— 120226 


SpeX Tft V 1 


— 120912 


SpeX 


— 121104 


SpeX 



A (yu.m) 



Fig. 14.— NIR/MIR spectra (2010 February to 2012 November) 
and photometry data (1996 March to 2008 March 26) revealing 
significant variability on time scales of several months to years. 



correlation between the equivalent width of the Br 7 line 
and NIR continuum infrared excess (Fig. [TBI bottom). 

Our multi-epoch spectro-photometric data from visual 
to MIR wavelengths are overplotted in Fig. [14] (and in 
Fig. EJfor our MIDI and BASS spectra), revealing signif- 
icant variability, which seems to have the strongest am- 
plitude (15%) in the K band, with a calibration uncer- 
tainty of 5%. Given that the amplitude of the continuum 
variability seems to peak in the K band and decreases 
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TABLE 3 
Results from our geometric model fits to the T-ReCS 

DATA 



Geometry 


A c 


FWHM 


-'"ring 

[mas] 


i 


PA 




[Mm] 


[mas] 


[°] 


[°] 


GAUSS 


8.8 


113 ±17 




40.8 ± 2 


119 ±9 


GAUSS 


9.7 


135 ±8 




31.8 ±7 


104 ± 15 


GAUSS 


11.7 


130 ±4 




31.5 ±6.5 


113 ±12 


RING 


8.8 




64 ±9 


38.2 ±2.5 


119 ±9 


RING 


9.7 




79 ±5 


31.3 ±7.5 


104 ± 15 


RING 


11.7 




75 ±1 


30.7 ±6.3 


114 ±13 



both towards shorter and longer wavelength, we consider 
it most likely that the origin of the variability is associ- 
ated with the inner disk, possibly with rotation-induced 
projection effects, as discussed in Sect. |U Photometric 
and/or spectroscopic variability h as also been found for 
other (pre-)transitional disks (e.g. iMuzerolle et al.ll2009l; 
Espaillat et ;all 120111: iSitko et all I2012t iFlahertv et all 
2012ft . 



3.3. Spatially resolved constraints on the gapped disk 
structure 

3.3.1. Geometry of the outer disk and gap region 

Our T-ReCS speckle observations resolve the disk 
structure at mid-infrared wavelengths of 8.7, 9.7, and 
11.7 /im. Providing a complete two-dimensional uv-plsaie 
coverage to baseline lengths of 5.5 m, these data are per- 
fectly suited to constrain the geometry of the large-scale 
disk and to determine the disk orientation and inclina- 
tion angle. A visual inspection of the six independent 
power spectra (Fig. [5]) reveals a significant elongation 
along PA~ 20° (corresponding to an orientation of the 
disk major axis along PA~ 110°). In the power spec- 
tra for the Si-3 filter, this elongation pattern is super- 
posed with another weaker signature, resulting in an ap- 
proximately cross-shaped power spectrum (Fig. [5] middle 
row). Given that this pattern appears only the in Si-3 
data set, which has the lowest SNR due to the atmo- 
spheric ozone absorption, we consider it likely that this 
signature is caused by instrumental artefacts, such as 
wind-induced vibration in the telescope. In the probed 
high-visibility regime (V 2 > 0.5), the visibility profile is 
rather insensitive to the detailed brightness distribution, 
which allows us to derive the disk orientation and incli- 
nation using a simple elliptical Gaussian model, without 
taking the inner disk structure into account. The mod- 
els include a photospheric flux component, although the 
contributions from this spatially unresolved component 
are practically neglible (< 2% based on our SED decom- 
position). 

The aim of our Gaussian model fitting is to quantify 
the elongation that is present in the 2-D power spectra 
(Fig. [5J . Using a Levenberg-Marquardt ^-minimization 
algorithm, we determine the best-fit full width half max- 
imum (FWHM), the position angle of the major axis of 
the elliptical Gaussian (PA), and the inclination angle 
(i). We performed the fit for each of our wavelength 
bands independently and we estimate the parameter un- 
certainties for each wavelength band from the scatter be- 
tween our two independent measurements. The deter- 
mined parameters are listed in Tabble[3]and consistently 
show an elongation along position angle PA=110° with 



an inclination angle i = 35° (i.e. closer to face-on). The 
Gaussian model provides a very simple and robust esti- 
mation for the elongation of the MIR emission compo- 
nent, but does not take the additional information pro- 
vided by the MIDI data into account, which reveals the 
presence of a compact spatial component. 

The choice for the second model geometry (elliptical 
ring + unresolved component) is motivated by our com- 
bined MIDI+T-ReCS data set, which shows that the vis- 
ibility function at wavelengths of 8.7 /xm drops signifi- 
cantly at short baselines (V 2 — 0.5 at < 5.5 m), but then 
stays rather constant at a level of V 2 « 0.1 from 46 m 
to 125m (Fig. [TH 3rd panel from top). It is important 
to note that this visibility profile is substantially differ- 
ent from the ones obser ved towards classic al T Tauri or 
Herbig Ae/Be stars (e.g. lKraus et alj 2008) and provides 
model-independent evidence for the presence of two spa- 
tial components, namely an extended component (which 
we represent with a ring geometry) and a compact com- 
ponent that is only marginally resolved even at hectomet- 
ric baseline lengths. Denoting the flux contributions of 
the ring component with -F r i n g, the flux of the unresolved 
component with F umcs , the total measured visibility will 
be given by Vtotai(-B) = (F r in g V Ting {B) + F unrcs ) / '(Fring + 
-Funrcs), where V r i ng (B) is the visibility profile of the ring 
component. At long baselines (where the extended com- 
ponent is overresolved, i.e. V ring ~ 0), Vtotai(-B) will 
converge towards F unY0S /(F r i ng + -Fu„ rcs ). According to 
our MIDI measurements, this convergence value is ~ 0.3, 
which leads us to fix F unrcs /(F r i ng + F unrcs ) = 0.3 in our 
model. The fractional width of the ring has been fixed to 
0.2. The best-fit parameters are listed in Tab. [3j where 
we denote the best-fit ring radius with i? r ing- 

This two-component model yields our best estimate 
for the disk inclination angle (i = 31.3 ± 7.5°) and disk 
position angle (PA=104± 15°). The existing Keck and 
VLTI long-baseline intcrferometric observations are not 
able to further constrain these parameters, but provide 
highly complementary constraints on the detailed object 
geometry, as will be discussed in the following sections. 

3.3.2. Geometry of the inner disk 

Our AMBER+ASTRA interferometric observations 
trace a compact emission component with a character- 
istic size of 0.90 ± 0.38 mas at 2/im (ring radius), as- 
suming a ring with fraction width of 20%, as adopted, 
for instance, bv lMonnier k Millan-Gabetll2002l This is 
more than 80 times smaller than the extended MIR com- 
ponent seen by our MIDI+T-ReCS MIR observations 
(79 mas at 12Lim). With measured squared visibilities 
down to V 2 = 0.6, the emission is clearly resolved both 
in the H and K bands. The visibility level drops slightly 
towards shorter wavelengths (Fig. [2j second panel from 
top). However, when plotted against spatial frequency 
(Fig. [151 2nd panel from top), the visibilities at all wave- 
lengths follow the same visibility function. This indicates 
that the visibility drop simply reflects the increase in 
angular resolution towards shorter wavelengths, but not 
a temperature gradient in the emitting structure, nor 
the presence of an optically thick inner gaseous emis- 
sion component, such as observed in the disks around 
some Herbig A e/Be sta rs (|Eisner et al . 2007; Kraus et al.l 
2008; Isclla et al. 2008). We also do not find any indica- 
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TABLE 4 
Results from the companion model fit to the 

NIRC2 DATA 



Spectral 
band 


Am 

[mag] 


P 

[mas] 


PA 

[°] 


H 
K> 
V 


5.73 ±0.31 
5.22 ±0.18 
6.05 ±0.35 


94.9 ±3.0 
42.8 ±5.4 
114.4 ±9.2 


61.9 ±1.9 
306.6 ±3.4 
328.6 ±4.4 



tions for an extended NIR scattered-light halo, such as 
observ ed for Herbig Ae/Be stars like R CrA (jKraus et al.1 
I2009al) . 

Unfortunately, the current uw-coverage of the AM- 
BER+ASTRA observations is not sufficient to measure 
the inclination angle of the inner disk independently. 
Therefore, we assume in our following modeling that the 
inner disk is coplanar to the outer disk. 

The characteristic size of the NIR-emitting region 
(ring radius 0.90 ± 0.38 mas or 0.34 ± 0.14AU at 
d = 385 ± 15 pc) is significantly larger than the ex- 
pected size of the dust sublimation radius (0.085 AU) 
as computed for black dust properties (cooling efficiency 
e = 1; sublimation temp e rature 1500 K) using Eq. 12 
in IDullemond fc Monnicr (2010), which includes back- 
warming effects. In Sect. 13.4.21 we suggest that this in- 
creased NIR size is due to contributions from optically 
thin material located at larger stellocentric radii. Alter- 
natively, it would be possible to reconcile the measured 
and expected sizes by assuming a population of smaller 
dust grains with cooling efficiency e ~ 0.2. 

The closure phases measured with AMBER are zero 
within the uncertainties of ~ 5° (Fig. [5]), which is consis- 
tent with a centro-symmetric brightness distribution on 
scales of p < 30 mas (r < 12 AU). 

3.3.3. Asymmetric structures in the gap region 

Measuring closure phases with sub-degree accuracy, 
our NIRC2 aperture masking observations are sensitive 
to asymmetries in the brightness distribution on angular 
scales well below the diffraction limit of the 10 m Keck 
telescope. The data reveal non-zero phases in all three 
wavelength bands, with varying significance levels. In 
order to explore the origin of these asymmetries, we fit 
a star+companion model to our data, where we treat 
the three wavelength filters separately. The error esti- 
mation and detection limits in our code have been care- 
fully tested and fine-tuned in various ear lier studies (e.g. 
iKraus et aIJl20lH: Kraus fc Irelan J20121) . The fitting re- 
sults are summarized in Tab. [4j while the corresponding 
best-fit phase plots and x 2 -maps are shown in Fig. [7] 

The K'-hsaid data provides convincing evidence for the 
detection of a companion at a separation of 0.043±0.005" 
from the star and towards PA=307 ± 3°. With more 
than 10cr, the detection is comparable to the significance 
level of other companion detections around transitional 
or pre-transitional disks with the aperture masking tech- 
nique. From the brightness ratio (AK = 5.22±0.18mag) 
and the 2MASS if-band magnitude (m K ^ = 7.41 ± 
0.03 mag), we estimate the absolute magnitude of this 
potential companion to Mk = 4.70 ± 0.29 mag for a dis- 
tance of 385 ± 15 pc. Based on the evolutionary mod- 
els from iBaraffe et all (|1998l ). we estimate the mass to 



M com p = 0.3 M Q at 6 Myr. 

Fitting our companion model to the H and L'-band 
data yields detections at ~ 3<r or more. However, for the 
following reasons, these detections are in conflict with 
the companion interpretation of the K'-bsaid detection: 

1. The companion positions in the H and L'-band 
best-fit model deviate significantly from the posi- 
tion derived from the i^'-band data (Tab.|4j. 

2. At the location of the X'-band detection, we ob- 
tain 3cr non-detections with magnitude limits of 
AiJii, n = 5.7 mag and AL[ im = 5.6 mag. These 
upper limits place strong constraint on the color 
of the putative companion, yielding an unrealistic 
red color of A(H — K) > 0.48 mag compared to the 
star+inner disk system. For a 0.3 Mq companion, 
we would expect to detect the companion with a 
color of A(H — K) — 0.26 mag. Therefore, in order 
to reconcile the companion scenario with the color 
constraints, we would have to assume that the color 
of the companion is reddened by inner disk shadow- 
ing and/or that the phase signal in the other bands 
is reduced by contamination due to disk emission. 

3. For the L'-band data, we find several solutions 
(Fig. bottom right), indicating that the com- 
panion model is not a good representation of the 
true object morphology. 

We interpret this as strong evidence that the measured 
asymmetries do not trace a single point-source, but dis- 
tinct asymmetric structures on different spatial scales. 
The model fits suggest that these asymmetric structures 
are located on spatial scales of ~ 40-110 mas (Tab.|4|) or 
~ 15-40 AU, i.e. within the gap region. This interpre- 
tation is in line with our disk model (Sect. |3T4.2[) . which 
predicts that the optically thin gap material contributes a 
significant fraction of the H-, K'-, and L'-band emission 
(Fig. 1161 top panel) . Any inhomogenities in the distribu- 
tion of this optically thin gap material introduces asym- 
metries that appear as a non-zero phase signal in aper- 
ture masking data, mimicking a companion signal. For 
complex disk structures, such as spiral arms, the ampli- 
tude and direction of the asymmetry changes with wave- 
length, as different wavelengths are sensitive to emission 
from different stellocentric radii. The non-zero phases 
measured by our aperture masking data likely trace such 
asymmetric structures in the disk gap. Given its lim- 
ited angular resolution, our aperture masking data is not 
able to constrain the detailed spatial structure of these 
asymmetric features, although this should become fea- 
sible with future long-baseline interferometric aperture 
synthesis imaging observations. 

3.4. Multi-wavelength interferometry + SED modeling 
of the disk structure 

In order to constrain the disk geometry quantitatively, 
we combine in this section our SED and NIR+MIR inter- 
ferometric data and fit it using a model of a geometrically 
thin dust disk. Our model assumes that the intensity I v 
at each surface clement in the disk can be described as 



B„(T(r))(l-e- r "), 



(1) 
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where B v is the Planck function and T(r) is the temper- 
ature at the considered disk radius, and r„ is the optical 
depth, which is related to the surface density S and the 
dust opacity k v at frequency v with 77, = k„£. The radial 
dust temperature is computed based on the thermal bal- 
ance between stellar heat ing and radiative cooling usin g 
the relation presented by iDullemond fc Monnierl (|2010f ) . 
This approach provides a first-order approximation of the 
intensity profile, sufficient for our purposes, but does not 
take the vertical disk str ucture, nor shadowing effects 
(e.g. Espaillat et al.l [2010T ) into account. These effects 
should be constrained in future studies using a more so- 
phisticated multi-dimensional radiative transfer and hy- 
drodynamics simulation approach. Assuming a point- 
symmetric brightness distribution, the model is also not 
able to reproduce the asymmetries that we have detected 
with our aperture masking observations (Sect. [373. 3j) . For 
the dust properties we u se opacities for silicat e and car- 
bon (graphite) grains bv iDraine fc Led (J1984D and mix- 
tures of these gr ain species. We ado pt the grain size 
distribution from iMathis et al.l (|1977fl , but tested also 
grains with a single grain size. 

Given our earlier, model-independent indications for 
an inner and outer disk component (Sect. 13.3. 2|) . we di- 
vide the disk into three regions, namely the inner disk 
(ID, extending from -Rid, in < r < -Rid, out)) the gap re- 
gion (i?iD.out < r < -Rod, in), and the outer disk (OD, 
extending from -R D,in < r < Rob, out), where R m ,in, 
-Rid, out, -Rod, in, and -Rod, out are treated as free parame- 
ters. For the inner and outer disk, we fix the surface den- 
sity to an (arbitrary) very high value, ensuring that the 
emission is optically thick. The surface density in the gap 
region S gap is treated as a free parameter with a flat sur- 
face density profile (£(r) = const.), as motivated by hy- 
drod ynamic simulations of gaps i n protoplanetary disks 
(e.g. lKlevllT99fll : IZhu et al.ll20Tl . The model is fitted 
simultaneously to the SED and the interferometric data 
and the five aforementioned parameters are fitted using a 
least-square minimization process. In order to avoid lo- 
cal minima, we vary the initial parameters systematically 
on a grid. Given the heterogeneous sampling of the SED, 
we binned the spectroscopic and photometric data on an 
equidistant grid in logarithmic wavelength space. In the 
fitting process, the VLTI/ AMBER, KI/ ASTRA, Keck- 
II/NIRC2, VLT/NIRC2, VLTI/MIDI, and Gemini/T- 
ReCS data are attributed the same weight. 

In addition, our model includes the stellar photosphere, 
which we represent with the recommended S YNTHE 
mode l for a FOV star (log 5 = 4.5, T cff = 7250 K,|Kurucl 
119931) We compute the flux ratio between the photo- 
spheric and circumstellar emission for each wavelength 
channel separately, yielding values from 45% at 1.5 /im, 
21% at 1.9 //m, to 12% at 2.5 /zm. 

3.4.1. Inner & outer disk (ID+OD model) 

A visual inspection of the SED (Fig. [T]) suggests that 
two simple blackbody components might be sufficient to 
fit the full SED. Therefore, in a first modeling attempt, 
we include only two optically thick emission components, 
without emission from inside the gap (S„ ap = 0). A 
standard dust mixture of 50% silicate + 50% carbon was 
used for this model. 

The resulting best-fit is shown in Fig. [15] and can re- 
produce the SED at NIR and far-infrared wavelengths 



Total model flux 
Inner disk 
Outer disk 






1 11.7(1111 . 


~!~-l 








■ ■ 1 ■ flW^ 



0.8 



0.6 



0.4 



0.2 



12 3 4 5 

Spatial Frequency [1 6 X] 

Fig. 15. — SED (top panel) and squared visibility am- 
plitudes (2nd to 4th panel), compared with the inner & 
outer disk model ("ID+OD" model, Sect. 13.4.11) . The 

data points at NIR wavelengths are color-coded as follows: 
blue: VLTI/AMBER+Keck/NIRC2, H band; red: KI/V2-SPR, 
K band; dark-red: VLTI/AMBER, upper K band; green: 
VLTI/AMBER+Keck/NIRC2, lower K band. The 2nd panel cov- 
ers the wavelength range 1.4-2.6 /im (H and K band) and the shown 
two curves cover the extreme values of this wavelength range. 
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reasonably well, even though it underestimates the flux 
in the MIR regime. The corresponding parameters are 
listed in Tab. EJ 

We find that this model results in much too high MIR 
visibilities, indicating that the MIR-emitting structure 
in the model is too compact. Improving the fit cannot 
be achieved by extending the inner disk (i?iD.in), since 
this would result in a conflict with the NIR visibilities. 
Increasing -Rid, out, on the other hand, changes the shape 
of the NIR excess and leads to an inconsistent MIR SED, 
suggesting that the optically thick inner disk is very nar- 
row. 

The main inadequacy of the ID+OD model is its fail- 
ure to reproduce the measured low MIR visibilities, in- 
dicating that this model lacks an emission component 
on spatial scales at least one order of magnitude larger 
than the inner disk. Therefore, we investigated whether 
parts of the outer disk might be able to contribute ex- 
tended MIR emission, for instance from a vertically ex- 
tended wall located at the inner truncation radius of 
the outer disk. In order to reproduce the MIR visibili- 
ties, this emission has to originate from stellocentric radii 
r > 46 AU, while temperatures T > 400 K are required in 
order to contribute significant amounts of thermal emis- 
sion in the 8 to 13 /jm wavelength regime. However, we 
were not able to find a physical scenario that would al- 
low us to achieve such high temperatures at the required 
stellocentric radii of r > 46 AU. For black dust prop- 
erties (e = 1), the expected temperatures at 46 AU are 
T ~ 80 K. This conflict cannot be resolved by adopting 
small grain species, such as 0.1 /im-sized carbon grains 
(e = 0.12), where the temperature would still be only 
T ~ 130 K. Therefore, we conclude that the discrepancies 
of the ID+OD model cannot be improved by changing 
the structure of the outer disk. This conclusion is also 
consiste nt with the finding from earlier SED-based stud- 
ies (e.g. iD'Alessio et~aLll2005t lEspaillat et~aTll2010D that 
investigated the detailed vertical structure of the outer 
disk wall of stars in Taurus and Ophiuchus and found 
that its contribution shortwards of 20 /im are negligle. 

3.4.2. Inner & outer disk with optically thin gap material 
(ID+GapMaterial+OD model) 

As a second scenario, we consider that the gap might 
not be completely depleted of material, but filled with 
optically thin dust. For this purpose, we vary both the 
surface density S gap and the dust properties. 

Our best-fit m odel is shown in Fig. [16] and includes car- 
bon dust grains dDraine fc LeeHl984Tl with a MRN-type 
grain size distribution ([Mathis et al.l I1977D . The dust 
mineralogy is constrained both by the absence of a 10 /im 
silicate feature and the shape of the SED, in particular 
the very sharp flux minimum around A ~ 15 /im (Fig.[T]). 
We find that small silicate grains are not compatible with 
our observations, since they result in a strong silicate fea- 
ture. Large silicate grains, on the other hand, result in 
too much excess emission at wavelengths A > 8 /im. A 
carbon dust composition with an MRN grain size distri- 
bution allows us to achieve a good fit with a \r total = 2.8 
(Table EJ. 

Compared to the ID+OD model, the NIR emission 
contains here contributions both from the inner disk and 
from the optically thin dust in the gap region, which re- 
sults in lower visibilities and, at the same time, a more 
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Fig. 16. — Same as Figure [15] but for the model with optically 
thin mater ial located in the gap ("ID+GapMaterial+OD" model; 
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TABLE 5 
Best-fit parameters for models described in Sect. 13.41 



Parameter 



ID+OD ID+GapMaterial+OD 



Inner Disk, inner radius 


-Rid, in 


AU 


0.19 


0.18 


Inner Disk, outer radius 


^?ID,out 


AU 


0.34 


0.27 


Outer Disk, inner radius 


^?OD,in 


AU 


44 


16 


Outer Disk, outer radius 


^?OD,out 


AU 


85 


85 


Dust surface density in gap 


S gap [g 


cm" 


- 2 ] 


9 x 10~ 6 



X r v 

V 2 
A r ,SED 

Xr, total 



73.5 
3.8 

65.6 



3.0 

1.7 
2.8 



Note. — X r v denotes the reduced x likelihood estimator for the visibility data, as 
defined in Kraus ct al. (2009a, equation 1). \ r sed i s defined equivalent for the SED data 
and x r total f° r t ne fit to the combined visibility and SED datasct. 




20071). or grain gro wth (|Dullemond fc Domimkl 120051: 
Birnst iel et al.l f2012f ). On the other hand, the observed 



Fig. 17. — Illustration of the V1247 0ri system, as constrained 
by our spatially resolved observations (sketch is not to scale). 

shallow visibility profile, providing a significantly better 
fit to the measured H and K-haxid visibilities (second 
panel from top in Fig. 1151 and [Tf5 j) . The improvement is 
even more significant at MIR wavelengths, where the N- 
band flux in the ID+OD model is dominated by contri- 
butions from the inner and outer disk (resulting in much 
too compact emitting structures at 8.7/im and 11.7/zm; 
third and fourth panel from top in Fig. I15|) . while in the 
ID+GapMaterial+OD model, they are dominated by op- 
tically thin material located in the gap region, resulting 
in a much better respresentation of the MIDI and T- 
ReCS visibilities (Fig.QU). 

4. DISCUSSION 

Our study provides spatially resolved constraints 
on the structure and physical conditions in a pre- 
transitional disk (Fig. [T7|) . enabling new insights on the 
disk clearing mechanism in these objects. We localize a 
narrow ring of hot material at the expected location of 
the dust sublimation radius, which rules out mechanisms 
that clear the disk in an inside-to-outside fashion, suc h 
as photoevaporation (e.g. lAlexander fc Ar mitagc 2007), 
magnetorotational instabilities (Ch iang fe Murrav-Clavi 



disk structure supports theories for dynamical disk clear- 
ing by planetary- or sub-stellar-mass companions, which 
might be responsible both for the gap-opening and for in- 
ducing the detected disk asymmetries. Our NIRC2 and 
AMBER observations rule out the presence of stellar- 
mass companions on scales from ~ 1 mas to a few arc- 
seconds, with brightness ratios up to 1 : 100 (correspond- 
ing to companion masses of 0.3 M Q ), but cannot rule out 
a very close spectroscopic companion. However, such a 
close stellar-mass companion would like also distort the 
inner disk, either by truncation or by causing asymme- 
tries, that have not been observed by our AMBER data. 

One of the main results from our study is the detec- 
tion of significant amounts of optically thin material lo- 
cated within the disk gap. SED-based models of several 
transitional disks also included an optically thin com- 
ponent, although in those cases the component has been 
introduced primarily to reproduce spe ctral features, such 
as the sil icate feature (e.g. TW Hya, ICalvet et al.l 120021 : 
LkCal5. lEspaillat et al.ll2011l ). Such spectral features 
are not present in our V12470ri spectra and we require 
the optically thin component in order to reproduce the 
MIR continuum visibilities. We find that this component 
even dominates at MIR wavelengths , which is in contrast 
to HD 100546 (iBenistv et all [20151 iTatulli et all l20ll 
and TWHya dRatzka et all 120071 : lAkeson et al.l I2011D . 
the other transitional/pre-transitional disks that have 
been extensively studied with infrared interferometry be- 
fore. No optically thin gap material was found in these 
objects, possibly indicating that TWHya and HD 100546 
are already in an later stage of gap clearing and/or grain 
growth than yi247 Qri. Adopting the dust opacities by 
iDraine fc Led (|1984h . we could achieve no satisfactory fit 
with silicate dust grains, but require a carbon-dominated 
dust mineralogy for the optically thin material located in 
the disk gap region. Of course, this conclusion depends 
to some degree on the adopted dust properties and future 
studies might test also other grain species. 

It remains an open question, how this carbon- 
overabundance has developed, but it is interesting to 
note that studies on the (3 Pictoris planetary system 
have also reveal ed a carbon-rich c omposition of the 
debris material (jRoberge et al.l 120061 ). Also, carbon- 
rich plariete_liave been predicted by theoretical stud- 
ies (jSeager et al. 2007) and indications for a carbon- 
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rich interior have already been found f or the super-earth 
55Cancrie ([Madhusudhan et al.1120121) . 

Using Keck/NIRC2 aperture masking observations, we 
measured strong non-zero phase signals, indicating small- 
scale asymmetries in the brightness distribution. Simi- 
lar phase signals have been observed on several other 
transitional and pre-tra nsitional d isks (e.g. CoKu Tau/4 : 
llreland fc Krausf l2008t T Cha: IHuelamo et all I201T1: 



LkCa l5: | Kraus fc Irelandll2012l; H D142527: iBiller et all 
120121: FL Cha: iCieza et all 120131 ) and have often been 
interpreted as signatures of a close companion. Fit- 
ting the -ftf'-band data, we also obtain evidence for a 
close companion (IOct detection significance; Fig. [71 mid- 
dle column). However, the direction and amplitude of 
the asymmetries changes with wavelength, which sug- 
gests that they are not caused by a close companion. 
We also rule out that the asymmetries might be caused 
by the heated wall of an inc l ined c entro-symmetric disk, 
as proposed by ICieza et "all (|2013|) , as such asymmetries 
would be consistently directed along the minor axis of 
the disk. Instead, our observations reveal more complex, 
radially extended structures, whose detailed geometry is 
not sufficiently resolved by our aperture masking obser- 
vations. We speculate that the asymmetries might be 
caused by spiral-arm features, similar to the ones discov- 
ered with Subaru/HiCIAO coronagraphic observations 
in th e pre-transitional disks of HD 1 35344B flMuto et al.1 
120121) and HD 36112 f=MWC 758. IGradv et al.l 12013ft 
on ~ 5x larger angular scales. In both cases, 
the asymmetries might resemble spiral density waves 
that are induced by orbiting sub-s tellar mass com- 
panions (jGoldreich fc Tremaind |1979() , although alter- 
native scenarios, su ch as gravitational instabilities (e.g. 
iDurisen et a l.l I2007T) or accr etion streams onto forming 
planets (e.g. iZhu et alJ l2012) have also been proposed. 

A major open question concerns the physical cause of 
the deep UX Ori-like occul tation events that were re- 
ported by ICaballerol ([2010). The detected two main 
events were separated by about 14 months and lasted for 
about two weeks each. Outsides these isolated events, the 
photometry of V1247 Ori was very stable on a level be- 
low 0.02 mag, which suggests that the occultations were 
not related to photospheric activity. ICaballerol ([2010D 
proposed that the photometric events might be caused 
by disk material passing the line-of-sight, assuming that 
the disk is oriented close to edge-on. At first sight, this 
scenario seems in contradiction to our inclination mea- 
surement of i = 31 ±8° (Sect. l3~3~l~T) . suggesting that the 
disks is seen closer to face-on than edge-on. However, 
it should be noted that this estimate has been obtained 
at MIR wavelengths for the outer disk regions, while the 
wu-coverage of our NIR data is still insufficient to ob- 
tain an independent estimate for the inclination angle 
of the inner disk. Therefore, we are not able to rule 
out that the inner disk might be considerably inclined 
or warped with respect to the outer disk. Hypothetical 
planets o n highly inclined orbi ts in the inner astronom- 
ical unit (jAlbrecht et al.l [20121 ) might shepard, deform, 
and tilt the inner disk with respect to the outer disk 
(jMuzerolle et al.1 [2009) , which might result in the ob- 
served occultation events. Rotation-induced projection 
affects of this inclined disk might also be responsible for 
the detected continuum variability, which peaks at NIR 
wavelengths and appears on time scales of several months 



to years (Sect. I3.2p . matching the orbital period of the 
inner disk. Alternatively, the occultation events might 
be related to the asymmetric structures that we have 
detected on scales of ~ 15-40 AU, for instance through 
a large, vertically extended disk warp that intersects 
the line-of-sight to the stellar photosphere sporadously. 
However, both for time scale arguments and hydrody- 
namic arguments, we consider the aforementioned sce- 
nario, which relates the occultation events with the inner 
disk, more likely. 

5. CONCLUSIONS 

We have presented coordinated near- and mid-infrared 
interferometric observations of a star with a pre- 
transitional disk, which has allowed us to investigate the 
disk structure of one of these intriguing objects in un- 
precedented detail. 

Our observations have targeted the previously sparsely 
studied pre-transitional disk around VI 247 Ori and re- 
solved a narrow, optically thick inner disk component 
located at 0.18 AU from the star, which dominates the 
near-infrared emission. This hot inner disk contributes 
less than 20% to the mid-infrared (8-13 /im) regime, 
while the largest fraction arises from optically thin ma- 
terial located in the gapped disk region from ~ 0.2 
to 46 AU. This strong mid-infrared continuum emission 
from carbonous dust grains located inside the gap re- 
gion has not been predicted by the current class of SED- 
based disk models. Multi- wavelength interferometric ob- 
servations on a larger sample of transitional and pre- 
transitional disks will be necessary to investigate how 
common this optical thin emission component occurs and 
whether it might be indicative of a particularly early 
stage of disk clearing. The organic bands (hydrocarbon 
features) originate from similar stellocentric radii as this 
optically thin material. 

Our aperture masking observations reveal highly sig- 
nificant asymmetries in the H-, K'-, and L'-band. The 
direction and amplitude of the asymmetries changes 
with wavelength, which leads us to reject a compan- 
ion origin and to exclude that the asymmetries are in- 
du ced by a hea t ed dis k wall, as proposed in the scenario 
by ICieza et "all (|2013l) . Instead, we propose that the 
asymmetries are related to strong density inhomogeni- 
ties in the gap region, possibly caused by the dynami- 
cal interaction of the gap material with the gap-opening 
body /bodies. Our results illustrate the highly dynami- 
cal and complex nature encountered in the inner gap re- 
gions and demonstrate the need for further high-angular 
resolution imaging studies. We would like to emphasize 
that the observed disk asymmetries can be easily misin- 
terpreted as companion signals and multi-epoch, multi- 
wavelength studies are essential for ruling out false de- 
tections. 

With a mass of 1.9 Mq, V1247 Ori is considerably more 
massive than other pre-transitional disks, approaching 
the intermediate-mass regime. Studying planet forma- 
tion in such intermediate-mass systems (spectral type 
A or B) is particular interesting, as abundance studies 
suggest that the planet forma tion efficiency increa ses to- 
wards higher stellar masses ([Johnson et al.1 12007). Be- 
sides being more abundant than those around solar-type 
stars, giant planets around early-type stars appear to fol- 
low also a different major-axis distribution, with a sig- 
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niheant deficit in the inner 1-2 AU, likely because of the 
different disk environments in which they formed. Multi- 
wavelength interferometric observations on V1247 Ori 
and other transitional disks now provide the exciting op- 
portunity to study the very early evolutionary phases of 
such systems and to explore the impact of planet forma- 
tion on the disk environment directly. Such observational 
evidence is essential to test planet formation theories and 
to link our knowledge of planet formation and disk evo- 
lution to the planetary system demographics observed in 
main-sequence systems. 
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